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ABSTRACT 

We describe the diffusion and random velocities of sohd particles due to stochastic forcing by 
turbulent gas. We include the orbital dynamics of Keplerian disks, both in-plane epicycles and vertical 
oscillations. We obtain a new result for the diffusion of solids. The Schmidt number (ratio of gas 
to particle diffusivity) is Sc w 1 + {Qtsto-pY, in terms of the particle stopping time tgtop and the 
orbital frequency i7. The standard result, Sc = 1 + tstop/^oddy, in terms of the eddy turnover time, 
^oddyi is shown to be incorrect. The main difference is that Sc rises quadratically, not linearly, with 
stopping time. Consequently, particles larger than 10 cm in protoplanetary disks will suffer less 
radial diffusion and will settle closer to the midplane. Such a layer of boulders would be more prone to 
gravitational collapse. Our predictions of RMS speeds, vertical scale height and diffusion coefficients 
will help interpret numerical simulations. We confirm previous results for the vertical stirring of 
particles (scale heights and random velocities), and add a correction for arbitrary ratios of eddy to 
orbital times. The particle layer becomes thinner for ioddy > 1/^^ with the strength of turbulent 
diffusion held fixed. We use two analytic techniques - the Hinze-Tchen formalism and the Fokker- 
Planck equation with velocity diffusion - with identical results when the regimes of validity overlap. 
We include simple physical arguments for the scaling of our results. 
Subject headings: Disks, Planetary Formation, Solar Nebula 



1. INTRODUCTION 

The stirring of particles by turbulent gas is a funda- 
mental problem in ffuid dynamics. It is especially impor- 
tant in studying the formation of planetesimals because 
the (non-Keplerian) motion of particles below about a 
kilometer in radius is strongly influenced by gas drag. 
Our goal here is to include, for the first time, the effect 
of Keplerian orbital dynamics on the random speeds and 
diffusion of particles. 

Circumstellar disks are thought to be turbulent. Tur- 
bulence drives accretion onto protostars by dissipat- 
ing Keplerian shear and transferring angular momentum 
radially outwar ds. The magneto-rotational instability 
(Balbus and Hawley 1991) is a leading candidate for the 
source of disk turbulence despite concerns that proto- 
planetary disks may be insuffi ciently ionized in "dead 
zones" near the disk midpla ne ( Gam mie '1996'; 'Semenov' 



et al. 2004[). Convection (Lin and Papaloizou 1980 



determines the Toomre ( 1964 ) gravitational stability, or 
Q, parameter for thin disks. By analogy with single com- 
ponent disks, it is argued that gravitational collapse re- 
quires densities above a critical threshold, and/or Q val- 
ues below unity. This overlooks the important role of 
gas drag in the collapse of solids, which remo ves such 
thresholds ( |Ward|[l976l pOOO| |Youdin| [2005a|bl ) . These 
works show that collapse mto rmgs over many dynam- 
ical times is possible even for nominally gravitationally 
stable values of midplane density and velocity dispersion. 
Radial diffusion plays an important role in this type of 
dissipative collapse. This work establishes formulae for 
vertical density structure, velocity dispersion and radial 
diffusion of particles in turbulent gas disks. These re- 
sults can be used to construct more realistic models of 
gravitational collapse of solids into planetesimals. 

Turbulent stirring also det ermines the relati ve speeds 
with which particles collide (Volk et aZ.||1980| hereafter 



msms ( jDubruUe et a/.l2005| [Rincon et aZ"t2007 , Lithwick] 
20071 to generate turbulence are being actively mvesti- 



Klahr al\ \m^ and pur ely hydrod yn amical mecha- \^^^) ^^^^ can be used to help determine when parti- 

M r . — — . ^ — — ■ cles coagulate, fragment, or melastically dissipate energy 

( Dominik and Tielens|1997 ) . Since our work does not in- 
clude the spatial dynamics internal to eddies we do not 
specifically address relative velocities. However for large 
solids with long aerodynamic stopping times the rela- 
tive speeds are just the RMS random spee ds (added in 
quadra ture). Other works giv e refinement s ( Mark iewicz 
et a?.|1991| ) and analytical fits ( Ormel and Cuzzi,2007| to 
the VJMR relative velocities, and investigate the impli- 
cations for coagul ation, given various assumptions about 
particle sticking (Mizuno et al. 1988). VJMR does not 
give results for particle diffusion and does not include 
orbital dynamics. 

Recent years have seen a remarkable advance in numer- 
ical simulations of the dynamics of solid particles i n tur- 
bulent gas disks ( e.g. Johansen and Klahr 2005 : Fromang 
and Nelson||2005| [Carballido et aL||2005( ). The resuTteoT 



gated. For this work, we remain agnostic as to the source 
of disk turbulence, as well as to its detailed properties. 

While turbulence is beneficial for accretion onto stars, 
it is thought to hinder the growth of planetesimals in 
disks. The gr avitational collapse of smaller solids into 
planete simals fSafronov |1969 Goldreich and Ward|1973] 



Youdin and Shu 2002) is opposed by turbulence ( |WeT 
^enscLilling 1980). Turbulent diffusion balances the ver 



plane (jCuzzi et al. 1993| hereafter CDC, 
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this work can be used to make detailed comparisons to 
the motion of test particles in these simulations. Simula- 
tions have now begun to inclu de the dynamical feedback 
of particles on gas dyna mics (jJohansen et al. 2006 Jo- 
hansen and Youdin 2007). These works find that stream- 



ing instab ilities (Youdin and Goodman 2005 , Youdin and, 
Johansen||2007 1 lead to the clumping of particles, an ef- 
fect our turbulent stirring models do not include. Other 
particle concentration mechanisms, in cluding long-lived 
vortices (de la Fuente Marcos and Barge 2001), turbulent 
pressure enhancements (Johansen et a/. , . 200 6 ^, or low- 
enstr ophy r egions at the dissipation scale ( |Hogan and 
Cuzzi||2007l are neglected as well 



TEe sta ndard result for particle diffusion was given by 
CDC (see J 4.1.1 1. We improve on this result in two ways. 
First we evaluate the statistical growth of particle separa- 
tions, instead of making the (incorrect) assumption that 
diffusion is proportional to kinetic energy. Second, we in- 
clude orbital dynamics, both vertical oscillations and the 
epicyclic motion that couples radial and azimuthal forc- 
ing. Since the CDC result was correct for small, tightly- 
coupled solids, many previous results arc unaffected by 
the error. 

This work is organized as follows. We introduce the rel- 
evant timescales of the problem in f|2] Our main results 
are justified with simple order of magnitude estimates 
in S|3] We review particle stirring with the Hinze-Tchen 
formalism in §4.1[ b efore extending it to include vertical 



oscillations in §4.2|and epicyclic oscillations in ^4.3 We 



confirm our results with the Fokker-Planck equation in 
^ Se ction [6] addresses neglected effects suc h as colli- 
sions (j ]6.1| , the "crossing trajectori es" e ffect ([6.2), and 
particle feedback on gas dynamics {{ 6.3 1. We discuss the 
implications of our results in SJT] 

2. TIMESCALES, DIMENSIONLESS PARAMETERS, AND 
PARTICLE SIZES 

Particle stirring in turbulent gas disks is characterized 
by three timescales: tstop, ^cddy, and l/Q. The parti- 
cle stopping time, tgtopi is the exponential timescale for 
the decay of relative motion with gas via drag forces, 
e.g. equation d7|, and tstop is an increasing function of 
particle size (for a given material density) as detailed 
in appendix IX] The eddy time, teddy, is the correlation 
time of turbulent fluctuations (defined formally in equa- 
tion 111. For Kolmogorov turbulence, ioddy is also the 
turnover time of the largest-scale eddies. Smaller ed- 
dies have shorter turnover times, and contribute less to 
the turbulent kinetic energy. The orbital frequency, i7, 
is also the vertical and epicyclic oscillation frequency in 
Keplerian disks, where 1/J7 is also the radial shear time 
(modulo 3/2). 

We define three dimensionless parameters, of which 
two are independent. 



St 



"^^stop 
^stop/ ^cddy 



Ts/Tc 



(la) 
(lb) 
(Ic) 



to describe the parameter space of turbulent stirring, as 
shown schematically in Fig. [l] The dimensionless stop- 
ping time measures whether particles are tightly (rg <C 1) 
or weakly (ts ^ 1) coupled to unperturbed circular gas 
orbits. The Stokes number, St, also measures particle 



coupling, but to turbulent fluctuations. The dimension- 
less eddy time, Te, describes the effect of orbital shear 
on eddies. If eddies are destroyed in a shear time, then 
Te ^ 1 , and some local MRI simula. tions support this 
choice (Fromang and Papaloizou 2006 1.^ Orbital dynam- 
ics plays no role in turbulence with Te <C 1, which is un- 
sheared and would require an energy source other than 
orbital differential rotation, perhaps convection. Turbu- 
lent eddies with Tq ^ 1^ are extended azimuthally by 
orbital shear. MRI turbulence might well give rise to 
To ^ 1 eddies as indicated by the extended azimuthal 



structure seen in MRI simulations (Hawley et al. 19951 
and per haps particle trapping seen m global MRI simu- 
lations ( Fromang and Nelson||2005 1 . 



Given our ignorance of the detailed properties of tur- 
bulence in protoplanetary disks, we consider deviations 
from the favored value of To = 1 (in which case Tg = St). 
By keeping To general, we also gain a physical under- 
standing of its effect. Appendix |B] describes a simple 
model for turbulent velocities with arbitrary To , which in- 
cludes basic dynamical considerations of shearing disks 
and is useful in making numerical evaluations. We in- 
clude general expressions for our results so other realiza- 
tions of turbulence can be considered. 

Previous works on particle stirring which neglect or- 
bital dynamics effectively assume Ts,Te <C 1 (regions I 
and II in the lower le ft co rner of Fig. [ij. Our study 
of ve rtical motions in |4.2| and epicyclic oscillations in 
§4. 3| includes all regions of Fig. [1] The shaded region m 



Fig. [T] indicates where ioddy is the shortest timescale in 
the problem (regions II and III) . Here the Fokker-Planck 
equation can model turbulent stirring with the inclusion 
of velocity diffusion, as done here (^ for in-plane motion 
and by CFP for the vertical settling problem. 

The vertical settling timescale, igott, of tightly coupled 
particles, Tg ^ 1, at the terminal velocity for stellar grav- 
ity is tsett = ^/i^Ts). The vertical oscillations of loosely 
coupled particles damp in tgott = ^stop- 
be combined to give: 

^^sctt — Tg -I- 1/Tg 



These results can 



(2) 

The hatched triangle in Fig. fl] (regions V and VI) shows 
where the eddy times exceecTparticle settling times. 

2.1. Aside on the Schmidt Number 

This paper obtains a new result for the spatial diffusion 
of heavy particles. Dp. We compare this to the spatial 
diffusion of gaseous contaminants, Dg, which depends 
only on the time-integrated correlation of turbulent ve- 
locities (equation 11). Following CDC we call the ratio 
Sc = Dg/Dp, the Schmidt number, and we generally ex- 
pect Sc > 1 because particles should not diffuse faster 
than the gas. This notation is somewhat unfortunate as 



^ Note however that |Promang and Papaloizou| |2006) use a local 
Eulerian measure of the eddy time, which they argue should be 
similar to the more appropriate Lagrangian measure. 

^ Even though long-lived vortices are often invoked to trap par- 
ticles, turbulence with Tc 2> 1 does not necessarily concentrate 
particles within eddies. To trap particles, a vortex must live for 
many circulation times (the time for fluid to circulate within a vor- 
tex). For example, in fjBlwe introduce a simple model for Tc 2> 1 
turbulence that is composed of vortices that live for a single cir- 
culation time. Even though their circulation time is much longer 
than the orbital time, these vortices cannot concentrate particles 
in their lifetime. 
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it differs from the standard fluid dynamical convention. 
In fluid dynamics Schydro = ^/Dg, is the ratio of the 
viscous (momentum) diffusion to mass diffusion (of the 
fluid itself for the typical single fluid case). This ratio 
is astrophysically important since it determines the dif- 
fusion of gas (or s mall particle) con t aminants relative t o 
the accretion flow ( [Stevenson [l990l [Turner et alJ^m^. 
Most numerical s tudies of particle ditt'usion ( [Jonansen 



and Klahr 2005] [Carballido et a?.[ [2005) measure the 
er for particles as the ratio Scgim = vlD-o, 



Schmidt number for particles as the ratio Scgim = v/L-p, 
where v includes the total accretion stress, with Reynolds 
and Maxwell stress contributions, since the turbulence is 
MRI-driven. Thus Scgim combines the effects of parti- 
cle dynamics and ang ular momentum transport. Indeed, 
Johansen et al. (2006 ) found that imposed vertical mag- 
netic flux signiflcantly affected Scgim- This explains why 
there have been reports of Scsim both above and below 
unity for tightly coupled particles, while we flnd (with 
CDC) that Sc = 1 in this limit. To isolate the effects of 
particle dynamics, and to compare directly with CDC, 
we opt for the CDC convention of Sc. Comparison of 
our results with numerical simulations requires measure- 
ment of both the diffusion of heavy particles. Dp, and 
Lagrangian tracers of the fluid motion, _Dg. 



Larger particles 



a 1 




"^s - ^^stop 



Fig. 1. — Schematic plot of the parameter space for particle stir- 
ring with eddy time vs. particle stopping time, both normalized to 
the orbital frequency. Dotted lines indicate where these To and Ts 
parameters are unity. The straight line of slope one demarcates 
Stokes number of unity, where particles are marginally coupled to 
eddies. The shaded region indicates where a Fokker-Planck ap- 
proach is valid because the stochastic forcing time, teddy i is short. 
In the upper hatched region, teddy exceeds the vertical settling 
time. Roman numerals label regions where timescales are ordered 
as indicated in Table 



3. PHYSICAL SCALING ARGUMENTS 

Our main results can be understood with simple phys- 
ical arguments in limiting cases. We assume isotropic 
turbulence with a RMS speed of 5Vg. The turbulent dif- 
fusion coefficient for gas is Z?g ^ 6V^ teddy 



3.1. No Orbital Oscillations: Tg ^ l,rc <C 1 

We flrst ignore orbital effects by taking To <C 1 and 
Tg <C 1 (lower left quadrant of Fig. [ij, so the particle 
response depends only on St. For St ^ 1 particles are 
well coupled to eddies. Thus the particle RMS velocities, 
6Vp ~ 5Vg, and diffusion coefficient, Dp ~ Dg, are the 
same as for the turbulent gas. 

For St ^ 1 turbulent eddies wash over loosely coupled 
particles, and the kicks give a random walk in velocity. 
Each kick of duration teddy has an amplitude T^kick ~ 
^V'g toddy /^stop in a random direction. Particle motions 
equilibrate with the turbulent forcing in tstop • Thus after 
N — tstop /teddy kicks, the random velocity will saturate 
at: 

5Vp w ViVViiek = (5Vg/VSt for St > 1 . (3) 

By combining the large and s mall St limits we recover 
the standard result of equation ( [13[ ) . 

To estimate particle diffusion for St ^ 1, note it takes 
tstop to deviate the particle velocity by order unity. Thus 
the particle mean free path is £p 



diffusion is 



'p ' - 5Vp £p ^ U Vp tstop ~ O Vg t-oddy 

With increasing St, the longer mean free path cancels the 
lower velocities to main tain Dp = Dg for all St in agree- 



SV^t, 



iJV'pistop, and the 



D. 



ment with equation ( 14 1, which also ignores orbital oscil 
lations. Analysis of previous results which erroneously 
claim that Dp dec reases for large St in this regime can 
be found in §4.1.1[ 

3.2. Including Orbital Oscillations: Tg 3> 1 3> To 

To better understand orbital oscillations we consider 
Ts 1, so that particles decouple from gas orbits, but 
keep Te <C 1. Thus particles receive many small kicks 
during a vertical oscillation or an epicycle and St 3> 1. 
This region occupies the lower right quadrant of Fig. [l] 
The particle velocity response is still a random walk. 



which saturates at 6Vp '--^ SVg/vSt as above, unaffected 
by oscillations with freq uency fi. This is conflrmed in 
equations ( 20 1 and ( 35 1 , where the latter includes cor- 
rections for the ratio of radial to azimuthal epicyclic mo- 
tion. 

Spatial particle diffusion is interrupted by epicyclic os- 
cillations with a length scale £opi ~ SVp/f2. Vertical os- 
cillations are centered on the midplane, and thus do not 
freely diffuse over long times. For epicyclic motion, we 
estimate diffusion by noting that the guiding center shifts 
by £cpi every tgtop, because the velocity changes by SVp 
in this same interval. This gives a diffusion coefficient: 



D, 



^cpi 



t 



stop 



svl 

]7ts 



eddy 



(4) 



This can be conflrmed by considering the more incre- 
mental changes during each orbital oscillation. Ev- 
ery 1/fI a particle receives A^epi = I/to velocity kicks 



ick 



SVg/St to give a velocity change, SVc 



cpi 



of 14i 

Vkick \/ Aopi ~ '^V'gy^/Ts, every orbital period. Thus 
every orbital period gives a shift in the guiding center 
position of 5^opi ^ (^V^pi/I?. A random walk with steps 
of size ^€opi every orbital time gives a particle diffusion 
coefficient Dp ~ '^^cpi''^ ^ -^g/'^s^ ^ i"^ equation (W| . 

If we include the Tg ^ 1 limit where orbital effects are 
negligible, then our estimate of the Schmidt number is: 

sc.§s^i+.; (5) 
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To order unity, this agr ees with the precise result for ra- 
dial diflFusion, equation ( 37 ) , which holds for all Tg. Thus 



orbital oscillations reduce the diffusion of loosely coupled 
particles. 

3.3. Vertical Scale Height 
We obtain the particle scale height, Hp, by equating 
the diffusion time, H^/Dp, and the settling time (equa- 



tions ^ and (|5|). The result. 



H„ 



(6) 



holds for arbitrary Ts and agrees with our detailed re- 
sults in equation (21 1 and previous works. In the loose 



coupling limit. Hp can be obtained simply by noting that 
Hp ^ SVp/fi for vertical oscillations. CFP computed Hp 
in the loose coupling limit both analytically and with 
hydrodynamic simulations and found equation (|6]).'^ The 
simple link between equations (jsl) and (6| confirms that 
orbital oscillations weaken particle diffusion. 

4. FOURIER SOLUTION OF LANGEVIN EQUATIONS 
4.1. Hinze-Tchen Formalism: No Orbital Oscillations 



W e summarize t he classic derivation of Hinze ( 1959 1 
and Tchenj (19471 for turbulent transport of particles 
in the absen ce of orbital effects or other external forces 
(reviewed in Fan and Zhu 19981. A Langevin equation 



describes the forcing of particle velocities, Up, by a gas 
velocity, Ug, with stochastic fluctuations as 



dup 
~(it 



Ug{t) 



is 



t. 



(7) 

^stop ''stop 

We restrict our attention to a single dimension with no 
loss of generality since the motions in different dimen- 
sions are uncoupled. The drag force is linear in the rela- 
tive velocity, which applies in either the Epstein or Stokes 
regimes. For protoplanetary disks, this requires particle 
sizes below aturb in equation ( A4 1 . Spatial variations in 



the gas forcing are ignored, because the ti me to cross an 
eddy is longer than the eddy lifetime (see [6.2 1. 

We decompose Up and Ug into Fourier modes that vary 
as exp(— zwt). The Fourier amplitudes evolve by equa- 
tion ^ as: 



1 — ZW^stop 

The power spectrum of turbulent forcing, 

-^g) ^cddy 



Egico) 



(8) 



(9) 



is normalized to the expectation value of turbulent ve- 
locity squared. 



dujE„ 



(10) 



The frequency dependence in equation ^ corresponds 
to exponential decay of time correlations, as the Fourier 
transform shows, P{t) = exp(— |t|/ioddy)/(27r). This 

CFP also show that both the spatial and velocity distributions 
are Gaussian with the stated RMS centered on zero velocity at the 
midplanc. 



spectrum is Kolmogorov since P(w) cx w ^ for frequen- 
cies above 1/icddy, the integral scale. This equates to 
the familiar spatial spectra of P{k) (x k~^^^ since kvf. = 
constant, and Vk — t-^/k readily gives k-^/^dk 
Diffusion of tracer particles (i.e. gas) with position, 

Xg{t) — J^Ug{T)dT + Xg(0), is given by the statistically 
averaged growth rate of squared displacements over long 
times as 



\d{xl 



2 dt 
dr 



(ug(r)ug(0))dr (11a) 



dujEg{uj)e^-"^^^ (lib) 

= ^^g(0) = {ul)t,My (lie) 



where eq. (11d) equates the velocity autocorrelation 
function to the Fourier transform of the power spec- 
trum, and eq. (jTl]:) exploits that the power spectrum 
is even in lo (equivalent to the time invariance of the 
correlation function). The final result is the expected 
product of RMS velocity SVg — (Wg)^/^ and eddy length, 

^oddy — ^^g. ^eddy- 

The power spectrum of particle motions 



Epico) = Egiu:)\ 



1 



(W^stop)^ 



(12) 



readily produces the mean squared particle velocities 



dujEp^Lj) 



{ul)t 



eddy 



eddy 



+ t 



stop 



l + St 



(13) 



This agrees with our physical estimate in equation ([3| , 
and more importantly, with a detailed treatment that 
includes the spatial spectrum of eddies (VJMR). 

Spatial particle diffusion is computed as in equa- 
tion (111. Since the power spectrum in equations ([9| 
and (T2| is also even in w, we simply have 



Dr. 



TTEpiO) 



(Ug) teddy — Dg , 



(14) 



independent of stopping time. The Schmidt number is 
Sc = Dg/Dp = 1 in this regime, as explained in *|3]by 
the longer mean free path for larger tstop- 

4.1.1. Comparison with CDC 

The most often cited result for the Sch midt number 
in the astrophysical literature is from CDC (fC uzzi et al. 
1993[) . Their analysis begins exactly as above (j j4.ip : us- 
ing the Hinze-Tchen formalism with the same temporal 
power spectrum and ignoring orbital dynamics (see their 
appendix B). Instead of calculatin g di ffusion fro m ve- 
locity correlations as in equations (111 and (14), they 



assume that diffusion coefficients are proportional to the 
ratio of fluctuating kinetic energies. This gives, ^ 



SccDC = {ul)/{ul) = 1 + St, 



(15) 



® CDC further add a correction for the "crossing trajectories 
effect" of particles drifting through eddies with a speed A V to get 
their final result of: Sc^^f. = (1 + St)(l -|- Ay2/^^2^)i/2^ n^^a^ 
further increases Sc and does not explain the discrepancy. See ^6.2| 
for further discussion. 
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TABLE 1 

Summary of Results and Ordering of Timescales by Region in Fig. [T] 



Region 


Timescales 


n In ^ 


I 


Ts < To < 1 


1 


II 


To < Ts < 1 


1 


III 


To < 1 < Ts 


4/ts2 


IV 


1 < Tc < Ts 


4/ts2 


V 


1 < Ts < Tc 


4/ts2 


VI 


Ts < 1 < Ts^"' < Tc 


1 


VII 


Ts < 1 < Tc < Ts^""" 


1 



Hp/H^ Sup/Sug'^ Swp/Swg'^ 



\/ az/T^T e 

•Jo-zjrs 



1 



V ^5tc/2tI 

2/Tsd 
Id 

1 




Note. — Results for 
velocities from eqs. | |36[ l, \21 
velocity) assume isotropic radial 



artic le r adial d iffus ion {Dp,x), scale heigfit {Hp) and radial (i5np) and azimuthal {&Wp) RMS 
(ISSb), and \20\ in the appropriate limiting cases. The in-plane results (diffusion and radial 
il and azimuthal speeds. 



which is incorrect since it ignores the role of the particle 
mean free path. There is no disagreement for low Stokes 
numbers where Sc ~ 1. 

To better understand the CDC result for Sc, and why 
it appears in the geophysical literature as well, we de- 
fine diffusion coefficients for arbitrary time, Dp ^(t) 



X^^p,g(''")^*p,g(0))dT, which evaluate to (see also 
[Zbl|[l^ : 



Fan and 



{ul)t, 



eddy 



''stop 



4op[l-exp(-t/i,top)] 



'■eddy 



Dg(i)= (Ug)teddy 



eddy 

[1 - exp(-i/Wdy)]^ 
[1 - exp(-i/tcddy)] 



(16) 
(17) 



The short time limit, t <C tstopj^eddy, recovers the CDC 
result while the long time limit produces the diffusion 
result, as demonstrated graphically in Fig. [2] The short 
time limit is of interest in laboratory experiments and 
in some geophysical contexts for the dynamics internal 
to the largest eddies. Astrophysicists however are nor- 
mally interested in the evolution of protoplanetary disks 
over many eddy times. Furthermore, the early phase of 
turbulent mixing does not describe diffusion per se since 
s eparations grow linearly i n tim e, not as ^/i. 



Schrapler and Henning (20041 reproduce the Schmidt 
number result of CDC by alternate methods, specifically 
a mean field theory approach, and they also neglect or- 
bital dynamics. The source of the discrepancy with our 
result is more subtle in this case, but it comes down 
to their choice of a (second order) closure approxima- 
tion. A more accurate modeling of the correlation time 
with their approach should give agreement with our work 
(Schrapler, personal communication). In any event sim- 
ple physical arguments should remove any doubt that 
our result (actually the classic Hinzc-Tchen result) gives 
the correct particle diffusivity for the limiting case of no 
orbital dynamics. 



Another result for Schmidt number. 



(^1969 



SCSaf = (1 + St)^ 

hereafter S69), 



has been attributed to Safronov 

originally by CDC. We nnd no such claim in S69. CDC's 
attribution of this result is apparently due to an exten- 
sion of S69's theory, which itself contains errors. Specifi- 
cally, S69 gives 6Vp ~ SVg/{l + St), which is incorrect. If 
this is inserted in equation ( 15 1, it yields the aforemen- 



tioned expression for Scgaf. ihe apparent similarity to 
our result with orbital effects [equation ^ with Ts and 
St switched] is illusory. 



We are thus fully confident that turbulent particle dif- 
fusion is independent of St in the absence of orbital ef- 
fects, and for the simple turbulence models considered. 
This comparison is intended to give confidence in our 
methods and clarify a confusing field, not to denigrate 
these important works. 



II 



^stop^'^eddy 




eddy 



Fig. 2. — Schmidt number, the ratio of gas to particle mass trans- 
port coefficients, versus time in the absence of orbital dynamics or 
other external forces. Different black curves choose different values 
of St = tstop /teddy In the diffusion limit of long times, Sc = 1 
for all St {grey sold line). The jCuzzi et aZ.] ( |1993[ l result {grey dot- 
ted line segments) is valid for t small cornpared to the eddy and 
stopping times, but does not describe diffusion over long times. 



4.2. Vertical Oscillations 

We now modify the Langevin equation to include the 
vertical component of the Keplerian gravitational field: 



dwp 



Wg(t) - Wp 



(18) 



"-stop 



where 



— dzp/dt is the vertical particle velocity and 
Wg the vertical turbulent forcing, which obeys the power 
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spectrum of equation ([9]).^ Equation (18 1 describes a 
damped oscillator that is stochastically forced. The 
Fourier response of the particle velocity 



gives the squared random velocities as: 



(19) 



1+St(l+T2 



(20) 



which reduces to the Q — result of equation pSl) for 
eddies much faster than the orbital time. Fig. [5] plots 
the vertical RMS particle velocities from the square root 



behavior fol 



of equation (|20|) . A physical explanation for the large 

j]4.2.2 The limiting forms of equa- 
below) are presented in Table 



tion (20 1 (as well as eq. 121 
[ijfor each region in FigT^ 




cO 1 



cO 
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0.01, 100 
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Fig. 3. — Vertical RMS speeds of particles, relative to gas turbu- 
lence, vs. dimensionless particle stoppin g ti me for several choices of 
the eddy turnover time, from equation l |20| . Stirring is most effec- 
tive for Ts <C 1 and Tc ~ 1. For large Ts and / or smaller Tc , particles 
decouple from eddies. At large Tc rapid orbital oscillations cancel 
the kick from eddies, see j4.2.2[ 

The particle scale height. Hp, plotted in Fig. |4j is cal- 
culated using Zp = iWp/uj as: 



duj\z. 



1 



P I ~ 12 



Str^ 



1 + St 



(21) 
(22) 



using Dg^z — {Wg)tcddy We compare our result with 
the standard value derived by Dubrulle et al. ( 1995 ) for 

With (lo?) instead of (n?). 



Fig. 4. — Particle scale height relative to gas vs. stopping time for 
several values of the eddy time with turbulent diffusion at a level of 
az = Dg^z /{H^Q) = 10~^. The particle layer becomes thinner for 
increasing particle stopping time. The layer also becomes thinner 
for eddy times longer than the orbital time and the settling time, 
l/{f2Ts), due to the cancellation of kicks as in Fig. [s] 

tightly coupled particles and by CFP for loose coupling 



rr(std) 
P 



(23) 



where = Dg ^/ {H^Q) is a dimensionless parameter 
for turbulent diffusion,^ and iJg — c^jfl is the gas scale 
height. Equation ( 21 1 reduces to equation ( 23 1 when ^ — > 
1, i.e. for To <C Max(l, l/rg) as in regions 17II, III and 
VII of Fig. [1] A physical description of the large eddy 
time behavior follows. For To = 1 our result becomes: 



1 



i7r, V 1 + 2t,, 



(24) 



The maximum discrepancy of \/2 for Tg 3> 1 is not very 
significant. 

Since particles are confined to a vertical potential well, 
we cannot measure long term spatial diffusion as in equa- 



tion (111. We show in j |3.3| that the particle diffusivity re- 
quired to maintain Hp agrees with the physical estimate 
of Sc in equation ([s]) . 

4.2.1. Applications and Extensions 
The ratio between velocity dispersion and scale height, 



St 
1 + St ' 



(25) 



reduces to the Keplerian oscillation frequency for loose 
coupling to eddies, St ^ 1. One might have naively 
guessed that the criterion would be Tg 1, but see ^J3] 

^ Note that this is not the same Q-parameter used for diffusion 
of angular momentum, and that turbulence is subsonic for Oz/tb = 

(wD/a < 1. 
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and S 4.2.2 for explanations of the surprising behavior in 
regions II and V, respectively. For St <C 1, collision rates 
will fall below the usual disk estimate of !7 times the 
vertical optical depth, Tp by VSt, since the decrease in 
density (with larger Hp) is not matched by the increase 
in RMS speeds. Since particle relative speeds are reduced 
below RMS speeds by a factor VSi for St < 1 (VJMR) 
the collision rate is 



t, 



coll 



St 



1 + St 



(26) 



for equal size particles, assuming that vertical speeds 
are at least comparable to in-plane motions. Note that 
this result is independent of the strength of turbulence 
(though for different particle sizes turbulence must com- 
pete with differential aerodynamic drift). 

We also compare the vertical settling speed for Ts ^ 1 : 
'^sctt = l^zl^stop ~ ^TsH-p to RMS particle speeds: 



Wsctt 



(27) 



For Tg ^ 1, the above is a small quantity (except for th e 
extreme case of region VI, where Wdy > ^sott, see ^4.2.2 1, 
so that random motions exceed the ordered settling that 
would occur without turbulence. 

We have assumed that both the particle stopping time, 
and the turbulent spectrum are independent of height 
above the midplane. This restricts the range of valid- 
ity of our results to Hp < Hg so that gas density, pg, 
is roughly constant. In stratified models. Hp < Hg is 
imposed by the hypothesis that the quantity diffused by 
turbulence is not the absolute particle density, but the 
partic le concentration relative to gas® (DubruUe et al. 
1995[ their eq. 28) . We suggest 



H. 



-1/2 



< 1 



(28) 



as a simple way to enforce the restriction Hp < Hg. 

The approximation (that turbulence acts to smooth 
concentration gradients) does not hold for larger parti- 
cles, which move through the gas. In practice this dis- 
tinction is not important since Hp/ Hg « ^ ajTg, <C 1 for 
Tg 3> 1 and a < 1. Thus heavy particles always settle 
within Hg in any event. See ^for more on the diffusion 
of heavy particles. 

4.2.2. Physical Understanding of Vertical Forcing for Large 

In S|3]we provide physical explanations for the velocity 
dispersion, diffusion, and scale height of particles in the 
Te <C 1 limit, i.e. the lower half of Fig. [T] We now 
consider To 3> 1, which is considerably more complex, 
particularly in regions V and VI where toddy > ^sott- We 
denote the particle and gas RMS speeds as 5wp and Swg, 
respectively. 

We first note that region VII behaves just as region 
I since particles are tightly coupled to both orbits and 
eddies, and the eddy time is shorter than the particle 



* See also |Chamberlaiii and Hunten||l987[ l, p. 90, for the same 
approximation in planetary atmospheres. I'he assumption is that 
turbulent diffusion mimics molecular diffusion, which is d riven by 
concentration Kra, dients (and also thermal gradients, see [Landau] 
[and Lifshitz|1959| eq. 58.14). 



settling time. Thus 5wp ~ (5uv and H^ ^ yj Dg/{f2Ts) 



in agreement with equations (20 1 and (21), in the limits 
St < 1 and Str2 = Wdy/^sott < 1- 

Region IV has loose particle coupling to eddies and 
gas orbits, and long eddy times. As in regions II and III, 
the velocity grows by a random walk. However since 
^cddy ^ , the kick from an eddy is largely can- 
celled by the many vertical oscillations in an eddy time. 
Thus the kick is only received over the last orbital cycle, 
^kick ~ Swg/Ts. After iV = tstop/^eddy kicks the veloc- 
ity saturates with Swp ~ M{ick\/^ = Swg/ y/fsTc. Ver- 
tical oscillations with loos e orbit al coupling then gives 
Hp = 5wp/Q. Equations ( 20|21 1 agree with these esti- 
mates in the relevant St ^ 1, St = St Tg ^ 1 limit. As 
a consistency check, note that i?p/4ddy ^ ^ 1 

so particle oscillations easily fit within a single large eddy. 

Since eddy times exceed settling times in regions V 
and VI, Hp is not determined by the usual diffusion 
arguments. Instead long-lived eddies suspend particles 
at Hp ^ 5wg/{QTs), by force balance between gravity, 
fP'H p, and drag, Swg/tstopi in agreement with equa- 
tion ((fTl) for St < 1 and r ^/St = TsTq > 1 (regions 



V anoVI) . Hubbard and Blackman ( 2006 1 give similar 
arguments, but without the caveats on the ordering of 
timescales. Though "pinned" by slow eddies, particles 
will be stirred by faster eddies with t'^ 



eddy 



< t 



eddy i 



and a speed 6wg '--^ '^^g y ^cddy/^^ddy (the transforma- 
tion of the Kolmogorov relation ve cx (^^^ to the time 
domain). For < 1 (region VI), t^^^y ^ t^ott = 
{nTs)~^, and tight coupling to these slower eddies gives 



Swg/y/TsTc. For Tg > 1 (region V), 
^sctt = ^stop, and the relevant eddy speed is 
5wg^/St. However the fast oscillations give an 
effective kick duration of as in region IV, and 

Swp ~ Sw'g/Ts ^ Swgl ^TsTq. There is no accumulation of 
kicks since each eddy lasts a stopping time, toddy ^ ^stop- 



"-eddy 
6w'g r 



Equation (20 1 confirms, in the St = Ts/to ^ 1 and 
TgTo ^ 1 limits, that ~ (5wg/ ^t^t^ in regions V 
and VI as shown in Table |TJ 

These physical arguments justify our new results for 
the scaling of particle scale height and vertical RMS 
speeds in the large regime. 

4.3. Epicyclic Oscillations 

We now treat the turbulent forcing of epicyclic oscil- 
lations. For small eccentricities, we consider the locally 
Cartesian radial, Xp, and azimuthal, j/p, positions relative 
to a guiding center with Keplerian orbital frequency, f2. 
Hill's linearized equations of motion with gas drag read: 



Xp — lip 

2/p = Wp - {3/2)f2xp 

tstop 

Slup Vg{t) — Vp 



Up — 2Qvp 



(29a) 
(29b) 

(29c) 
(29d) 



stop 



where we measure azimuthal velocities, Up, relative to 
local Keplerian shear. We ignore any time-steady gas 
velocities, such as pressure-supported sub-Keplerian ro- 
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tation of the gas would provide.^ 

We assign the same temporal turbulent spectrum, 
P{u!) as equation For generality we will keep the am- 
plitudes of the radial, azimuthal, and correlated gas forc- 
ing; (Ug), ('^g), and (ugVg); independent in many expres- 
sions to allow comparison with numerical simulations, 
equation (Bl I gives a model to relate these quantities. 



4.3.1. Velocity Response 

The Fourier responses of particle velocities to turbulent 
forcing are: 

- - (1 - -s)% + 2r.*g ^3^^^ 



where Vs = wigtop- 
correlations read: 



(1-ZZ^s)2+t2 

-(rs/2)ug + (1 - iVs)vg 
The squared Fourier amplitudes and 



(30b) 



1 



d{uj) 



1 



At. 



1 

d(uj) 



UgV* + U*Wg 



UgV*g + iilvg 



1 



2tJv„ 



iLgVl + ulVg 



(31a) 



(31b) 



(31c) 



TiY - A{v^T^Y In eqns. (31 i-b) 



Mpl)* -I- MpUp 



where d{Lo) = [1 

terms in the numerator oc iVs{ugVg — UgVg) are dropped 
since they are odd in Vg and will vanish on integration 
over frequency. 

Since the amplitudes have the form, e.g. |upp = 
a\ugY + b\v gY + c{u gV*+ UgV g)/2, we generalize the par- 
ticle power spectrum from equation ( 12 1 to: 

2 



p^^iu;) = [a{ul) + b{vl) + c{ugVg)] P{uj) 



(32) 

and similarly for the y power and x~y correlations. 
The particle energies are derived as before, e.g. (Up) = 

I^oo ^P.xdi^, to give: 



^{[l + St(l + r,V2)] {ul) + 
2Tf{2 + St){vl)+2Ts{2 + St){ugVg)^ (33a) 
^1 r^(2 + St) + [1 + St(l + /2)] (.^) - 
r,(2 + St)(ugt.g)/2} (33b) 
{uD +T,i2 + St) {vl) + 
St-T^){ugVg)} (33c) 



[UpVp) - I ^ 



d2 
(I 



^ The resulting (locally constant) drift of particles relative to 
gas can s till be included with a correction to the eddy time from 
equation KOl . 



where d2 = {I + r^) [(1 + St)^ + r^] . 

Since the expressions are fairly cumbersome it is fruit- 
ful to investigate limiting cases. Table [T] lists the limit- 
ing forms of Sup/Sug = {{u'^ ) / {Ug)y^^ , using the turbu- 
lence model in equation ( |B1| for the other components of 
gas forcing. For more general analysis, we first consider 
Tg ^ 1, i.e. the left half of Fig. [T] The equations simplify 
considerably to: 



(Wp> _ _ (UpVp) 



(ul) {vD (ugVg) 1 + St 



(34) 



There is no epicyclic motion since particles are tightly 
coupled to gas orbits, and the response scales directly 
with the forcing. ^'^ The scaling with St agrees with equa- 
tion ( 13 1 in the absence of orbital dynamics. 



Next consider the limit of short eddy times. To ^ 1 and 
St 3> 1 in regions II and III, which we will compare to our 
Fokker-Planck results in f|5] Here the particle response 
is, 



(l + 4) {ul)+2T^{vl) + 2n{ugVg) 



St(l+T2) 



St 



1 + 5r2/2 



1 



(35a) 



,Vg)/2 



St(l + r2) 



1 + 5t2/8 



St 



■^{ul)+T,{vl) + (UgVg) 



St(l 
3r,/4 



St 1 



(35b) 



(35c) 



where the final equalities assume isotropic turbulence 



with 



and ignore 



consistent with 



equation (Bl ).^^ For small Tg equations (351 recover the 
St ^ 1 limit of equation ( 34 1 . For Tg ^ 1 we recover 
epicyclic motion. The radial RMS speeds are twice the 
azimuthal RMS speeds, whether the forcing is radial, az- 
imuthal or combined. ^2 Equation (35") shows that ra- 
dial turbulent forcing produces a negative (upVp), i.e. 
inward transport of particle angular momentum. Az- 
imuthal forcing has the opposite sign, and dominates for 
our assumption of isotropy. 

Region IV exhibits epicyclic oscillations since i7^^ is 
the most rapid timescale. The ratio of RMS particle to 
turbulent speeds scale at ■y/St/r g — l/ ^/rgTc in agreement 
with the physical arguments in §4.2.2[ Region V is rather 
complicated since even though Tg > 1 particles do not 
complete epicycles as they are tightly coupled to eddies, 
St < 1. 

A possible exception occurs if, e.g. the azimuthal turbulent 
forcing is much stronger than the radial. 

We have dropped a term — TaTn( up-v^) from the numerator of 
the middle expression in equation \35\ :) because it is always smaller 
than the other terms in the numerator for To <^ 1. 

Random velocities are measured relative to the local circular 
speed. Relative to a fixed guiding center the elongat ion of epicy- 
cles is opposite ; azimuthal speeds are twice radial ( [Binney and| 
|Tremaine |1987| . 
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Fig. 5. — Particle RMS speeds (radial, solid curves; azimuthal, 
dashed curves), relative to the turbulent RMS speed, Vturb = 

,J{vi^^~+v^^)j2, versus stopping time for three different values of 

the eddy time (teddy = 1-0/17, 0.1/17, and 10.0/17 in the top, mid- 
dle, and bottom panels). The (square root of) the particle velocity 
correlations is plotted {dotted curves). The turbulent gas correla- 
tions (grey lines) transport angular momentum, but have a small 
effect on the particle velocities (and only near = 1). See i ]4.3.2| 
for discussion. 



4.3.2. Numerical Evaluation and Visual Interpretation of 

Velocities 

Fig. 5] plots RMS velocities and velocity correlations 
from the (square root of ) eq uations ( 33 1 for the turbu- 
lence model of equation (Bl I. This nioael assumes that 
velocities are isotropic, 



was seen for vertical speeds in Fig. |3] Thus the RMS 
velocities of large particle are a good indicator of whether 
the eddy time differs from unity. 

The speeds near Tg « 1 are sensitive to the eddy time 
and the {ugVg) correlations. For Tq = 1.0, the "bump" 
in radial speed near Tg = 1 is due to the high value of 
(ugVg) (otherwise the radial velocities would be flat up 
to Tg — 1). For To = 0.1 there is a noticeable drop in 
speeds by Tg — 1, because with St = Tg/re = 10 parti- 
cles have begun to decouple from eddies. For Tq = 10 
the radial velocity bump exists even without gas corre- 
lations, which are weak here. The difference in speeds 
near Tg w 1 discriminates between large and small eddy 
times. We emphasize that the turbulent transport of an- 
gular momentum by Reynolds stresses {ugVg) has only a 
small effect on particle speeds, and then only near Tg « 1. 

4.3.3. Radial Particle Diffusion 
Radial diffusion ^'^ of particles can be calculated as in 



equation ( 14 1 to give 



dT 



eddy 



-4Tg^(^g) +4rg(Mg«g 

(l+Tg2)2 



(36) 



The Tg ^ 1 limit recovers the result of equation ( 14 1 that 



particles of any size diffuse as a gas tracer when orbital 
motions are neglected. Azimuthal forcing has a negligible 
effect since tightly coupled particles quickly adjust to the 
local orbital speed. For Tg ^ 1, radial diffusion is domi- 
nated by the azimuthal forcing, Dp^x — A{vg)tcddy / , 
since angular momentum perturbations are needed to 
change semimajor axes. The scaling for loose coupling 
agrees with the estimate of equation ^ . Correlated gas 
fluctuations, (ugVg), are a minor correction to radial par- 
ticle diffusion, and only near Tg ^ 1 (a s w ith the veloc- 
ities). The limiting forms of equation (|36|) are listed in 
Table [l] for the turbulence model in equation ( B 1 ) . 
The radial Schmidt number can be expressed for ho- 



mogenous isotropic turbulence 
as: 



SCx 



D„ 



(1 



D 



l+4Tg2 



= 



(37) 



{Ug) — (Vg) , and estimates a max- 
imum correlation, {ugVg) (grey curves). Speeds are nor- 
malized to the characteristic turbulent speed, Vturb = 
{{Ug + Vg)/2y/^, and results are shown for To = 0.1, 1.0, 
and 10.0 (middle, top, and bottom panels respectively). 

For Tg ^ 1, particle RMS speeds simply follow the gas 
forcing, here the isotropic value Vturb- For Tg ^ 1, the 
RMS speeds fall as Tg , and the ratio of radial (sohd 
curves) to azimuthal (dashed curves) speeds approaches 
the epicyclic value of two. As To deviates from unity 
(becoming either larger or smaller) the RMS speeds at 

large Tg drop roughly as Tc/^I + T^). Similar behavior non-axisymmetric structure 



where Dg^ — (Ug)ioddy- This is consistent with, but 
provides order unity corrections to the rough estimate 
in equation It is interesting that this result is com- 
pletely independent of the eddy time- 

Fig- l6] plots the radial particle diffusion, Dp x, rela- 
tive to the diffusion of gas, Dg^^- The black curve s us e 



equation (36 1 and the turbulence model in equation (Bl I 



Since the correlations (ugVg) are a minor contribution, 
the blac k curves overlap with each other and with equa- 
tion ( [37| (or rather its inverse, not plotted)- The overlap 
only occurs because (Ug) = (Vg) in all cases- If the veloc- 
ities were anisotropic then the diffusion for Tg 3> 1 would 
vary as (Ug)/(Ug), i-e- Dp^^ would be larger (smaller) if 
azimuthal kinetic energy exceeds (falls below) the radial 
contribution - 

We do not here consider azimuthal diffusion. It is less in- 
teresting since orbital shear is much more effective at spreading 
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The CDC value for Sc, equation ( 15 1 , is plotted for ref- 
erence (dotted grey curve). The physical shortcomings 
of this result (neglect of epicyclic osci llation s and par- 
ticle mean free path) are discussed in [4.1.1 While the 
CDC value over-predicts diffusion for very large particles 
Tg ^ 1, it tends to under-predict diffusion for marginal 
coupling, Ts ^ I. 
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at t = 300i7 ^, in agreement with the diffusion seen in 
the figure. 
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Fig. 7. — A Monte Carlo simulation of the radial diffusion of three 
particles undergoing epicyclic motion with Tb = 10^ and Tq = 10~^. 
The radial coordinate x is normalized to th e RM S turbulent gas 
speed divided by the Kepler frequency. See §4.3.4| for details. 



Fig. 6. — Radial particle diffusion, relative to gas diffusion, v ersu s 
stopping time. The black curves [from equations \36\ and JbTJ] 
overlap despite different choices of eddy time. For tight coupTmg 
{ts "C 1) particles diffuse as effectively as gas. For loose coupling 

(ts 3> 1) particle diffusion falls as t^'^ . The CDC result, Sc = 
1 -I- St with Te = 1 {grey dotted curve) overestimates the diffusion 
of large (ts ^ 1) particles and slightly underestimates diffusion of 
marginally coupled solids. 



4.3.4. Monte Carlo Diffusion Simulation 



Testing our results with detailed turbulence simula- 
tions (or experiments) is beyond the scope of this paper. 
Instead we visualize the combination of diffusion and 
epicyclic motion with the Monte C arlo simulation shown 
in Fig.l^ We integrated equations (29) numerically, with 
Ts = 1(P and Te = 10^^ (region IlTof Figure nj). As a 
crude model for the forcing by the gas, we set Ug(t) = ±1 
and Vg{t) = ±1 for a time approximately teddy More 
precisely, we updated Ug and Vg at every time step with 
probability d^/tcddyi where dt is the time step, and the 
decision to update was implemented with a random num- 
ber generator. If an update was warranted, the sign of Ug 
and of Vg was each randomly selected. The three curves 
in the figure show three realizations of the random num- 
ber generator. For the oscillation amplitudes, equation 



i35k) predicts that (u^) = {5/2){ul)Te/Ts = (0.016)^, 



implying that the peak-to-peak amplitude of x should 

in agreement with the oscil- 
lations seen in the figure. For the diffusion, equation (36 1 
predicts the diffusion constant -Dp,^ — tcddy{4:/Ts){v'^y^ 
r2-i(0.002)2, implying that x ~ {2Dp,jy^^ = n^^O.Ob 



5. FOKKER-PLANCK EQUATION 

When the eddy correlation time teddy is the smallest 
timescale, the random gas motions give short uncorre- 
lated kicks to the particles which can be described by 
a Fokker-Planck equation. The Fokker-Planck equation 
has been used exten sively to study Brownian motion 
( Chandrasekhar|1943 1 . While the uncorrelated collisions 
in Brownian motion lead to thermal equipartition, the 
correlated motion of turbulent eddies excite particles to 
much higher speeds. CFP use the Fokker-Planck tech- 
nique to study vertical settling. By solving the in-plane 
Fokker-Planck equation we reproduce the results of §4.3| 
in the small teddy limit, adding confidence in our results. 

It greatly simplifies the analysis to switch to indicial 
notation, re- writin g th e equations of motion for a particle 
given in equation ( 29 1 as 

3 

Ui - -HxiSiy (38a) 



dxj 
~dt 
dui 
~dt 



A, 



t 



stop 



(38b) 



where i,j G {x, y}, summation over repeated indices is 
implied, and 



A, 



stop 

J7/2 



2f2 

'•stop 



(39) 



We suppress the subscript p from Xi,Ui. 

The Fokker-Planck equation for the particle distribu- 
tion function f(xi,Ui,t) reads 
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2 UiUj 



(AwjAitj) 
At 



/ (40) 



(e.g., equation 8-53 in Binney and Tremaine||1987 l. The 
angled brackets denote an ensemble average over the 
stochastic variable Ug^i, and Aa;^, Am^ denote the change 
in Xi , Ui over the time interval At, in the limit that 
At 0. (To be more precise, At must be larger 
than toddy I but much smaller than all other timescales). 
We have dropped the term containing {AxiAxj)/At be- 
cause, from equation (38i), it is proportional to At, and 
hence vanishes for At 0. To calculate {AuiAuj) / At, 
we may drop the AijUj on the right-hand side of equa- 
tion (38 because it only contributes a term of order At 
to {AuiAuj) I At] we are left with 



Aui = 



stop 



UgA{t')dt' 



(41) 



implying 

(AuiAuj) 
At 



At t2 

'-stop 



t+At pt' 



2teddy 
t^ 

''stop 



{Ug4t')Ug,,{t"))dt"dt' 



(42) 



in agreement with eq. (9) of CF P. The remaining aver- 
ages are given by equations (|38|) as 



(Aa 



At 
{Auj 
At 



= (ui) - -n{xi)6. 



ly 



= A,. 



The Fokker-Planck equation then becomes 
suppressing angled brackets and defining 



u ■) 

''stop 



(43a) 
(43b) 

(44) 
(45) 



as the velocity diffusion, and dt = dt — {3/2)f2xdy which 
reduces to a simple time derivative for an axisymmctric 
{dyf = 0) distribution function. 



We solve the Fokker-Planck equation (44 1 by taking 
velocity moments, defining the mass per unit area, mean 
velocity, and peculiar velocity tensors: 



fdu^duy 



n,; 



Uifduxduy 



(uj - Ui){uj - Uj)fduxduy 



(46a) 
(46b) 
(46c) 



{ui - Ui){uj - Uj){uk - Uk)fduxduy{A6d) 



Multiplying equation (44) by {l,Ui,UiUj}, and then in- 
tegrating out Ux,Uy, and y yields 

dtS + diSUi = 0, (47a) 
dtU, + U.djU, = ~S-'d,{SIl,^) + AjU, , (47b) 

S-^^k{SEk^J) = A^feHfej + A.-fcHfe, + 2i?„.447c) 



We henceforth consider the axisymmctric case {dy = 0) 
for simplicity As shown below, the left-hand sides of 
equations (47 1 and (47 1 are typically negligible. Assum- 
ing for now that this is true, equation (47 1 can be solved 
for the velocity dispersion. 









— ^stop 1 


K^yyJ 





''Stop 



'l + T^/2 2t,_ 
-r,/4 1 




(48) 



which recovers 



tion(45l 



-r,/21 + r2/2, 
equations (|35| once we identify 
4) 



{{uD , {upv'^{vl)} and apply equa- 



Equation ( 47 1 can be solved for the mean velocity 

UyJ I^l + T^ \-rs/2 1 J ^\SUxy 



.(49) 



We now assume that S varies on a shorter length scale 
than tgtop or Ilij do. The diffusion equation, from equa- 



tion (47 1, takes the simple form 



where 



t, 



stop 



''Stop 



1 



(50) 



(51) 



in exact agreement with equation ( 36 1 . 

We now determine the conditions under which the left- 
hand sides of equations (|47b,c) may be neglected. Let L 
be the length scale (in the a;-direction) over which \xiS 
varies, and let T be the timescale. The n U j also varies on 
length- and time-scales^ L and T (eq. |49]), with ampli- 
tude Ui 
with n 



47 



), and Ylij is constant (eq. 48 ) 



L/T (eq. _ _ 
C^, where C is the typical random speed. 
The left-hand side of equation (47 1 is smaller than the 
first term on the right-hand side when L/T <C C. As 
long as this inequality holds, the two terms on the right- 
hand side tend to equalize, with C^T jl? ~ Aij, and the 
terms on the left-hand side of equation ( [47| (other than 
the Sfcij term) are smaller than the "All" terms on the 
right-hand side by (L/CT)^ <C 1; the Sfe^j is smaller than 
the "An" terms by L/CT <C 1. In sum, our approxima- 
tions are valid as long as L/T <C C, i.e., that the mean 
speed of the particles is less than their random speed. 
As long as this holds, the timescales for both the mean 
speed and the random sp eed to reach their quasi-steady- 
state values of equations (48 1 and (49 1 arc much shorter 
than the diffusion timescale T. 



We return briefly to the discussion in \ 4.2.1 of whether 
particle concentration relative to gas is the fu nda mental 
quantity diffused by turbulence. Equation (49) shows 
that, for Ts ^ 1 (to ignore the complexities of orbital 
motions), STixx is the diffused quantity, which is not 
proportional to relative concentration. However this re- 
sult is restricted to tstop ^ ^cddy The standard assump- 
tion that diffusion levels concentration gradients holds 
for small particles tightly coupled to eddies. Since this 
limit is not treated by the Fokker-Planck approach, there 
is no inconsistency. 
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6. NEGLECTED EFFECTS 

6.1. Collisions 

Interparticle collisions are neglected both to simplify 
the analysis and because they are often a small correc- 
tion to the dynamics, when compared to gas drag. The 
collision time is 

2p,a 

icoll = 7, (5^) 

where v^^i is the relative velocity between particles and 
yO, ~ 1 - 3 g/cm'^ is the internal solid density. For col- 
lisions between two species i and j/* the collision time 
for species i due to j is set by the space density of j and 
particle radius, a — > 2aj/{af + a|). 
If we compare to the stopping time in the Epstein 
p,a/(pgCg), we see that 

,Ep 



,En 

regime, tj^^ 



stop 



Pp Wrol 



< 1 



coll 



(53) 



where inequality follows since both fractions are less than 
unity in a smooth turbulent disk.^^ The particle density 
should be less th an t he gas density for our test parti- 
cle approach (see §631) and the relative motions between 
particles (whether induced by turbulence or aerodynamic 
drift) are always subsonic. 

In the Stokes regime (when the particle radius exceeds 
the gas mean free path. A, but the Reynolds number of 
the flow around the particle remains below unity) colli- 
sions become relatively more important as the stopping 
time increases by 4a/(9A) > 1. However turbulent drag 
sets in before collisions dominate in gas rich disks. For 



turbulent drag. 



iturb 
'•stop 



p,a/{pgAV), where AV is the 



relative velocity of particles and gas. By this point (Ap- 
pendix [a]) particles are sufficiently decoupled from the 
gas that Vi-ei ^ AV, and 



jturb 
"-stop 

^coU 



Pp 

Ps 



(54) 



gives the result (for turbulent drag) that drag is more 
significant than collisions if the local mass in gas exceeds 
that in particles. We showed above that this criterion 
is over-restrictive in the linear drag regimes (Stokes and 
Epstein). 

6.2. Spatial Spectrum of Eddies and Crossing 
Trajectories 

UnUke the detailed treatments pioneered by VJMR, 
we neglect the spatial spectrum of turbulence. This is 
valid if the time to cross an eddy, icross — ^oddy/^^ 
(with ^cddy the largest eddy scale and AV the relative 
motion of particles and gas) , is small compared to either 
the eddy time (so that a particle never sees the spatial 
extent of an eddy) or the stopping time (so that a particle 
is fully accelerated before crossing an eddy) . Thus spatial 
homogeneity is a good approximation for 

Across > Min(tstop, teddy) • (55) 

Our stirring model assume a single particle size. Since we 
work in the test particle limit with no particle collisions, our results 
can be applied to a dispersion of particle sizes with no additional 
approximation. 

If significant clumping to pp > PgCg/f^ei >> pg occurs, 
from e.g. streaming and/or gravitational instabilities, then colli- 
sions would become dynamically significant. 



If AV arises from turbulent forcing, then this criterion 
is satisfied. For tight coupling, St <C 1, tight coupling to 
eddies gives AV <C 5Vg. Thus tcross/^cddy — SV^/AV ^ 
1, and tcross/istop is even larger, so equation (55 1 is read- 
ily satisfied. For loose coupling, St ^ 1, me particle 
RMS motions are small so AV ~ SVg, the gas forc- 
ing. Thus tcross/icddy ~ 4ddy / ( i^^g Wdy ) = 1, and equa- 



tion (551 is only marginally satisfied 

Orbital oscillations further help to confine particles to 
an eddy because the size of epicyclic or vertical oscilla- 
tions, £cpi, is smaller than the eddy scale: 



<5K 



<5K 



1 



^cddy 



> 1. 



(56) 



toddy SVg VSt 

The above assumes Tg 1 and St 1 to that particle 
complete epicycles. We conclude that spatial homogene- 
ity is a safe approximation for pure turbulent forcing, 
with (as usual) possible corrections near Ts = 1. 

However if AV is dominated by net motion of particles 
relative to gas eddies, the spatial extent of the eddies can 
be significant. CDC follow Csanady (1963) in describing 
this as the "crossing trajectories effect" (hereafter CTE). 
Azimuthal and radial drift of particles due to radial pres- 
sure gradients, dP/dr, is the most likely source of net 
motion. The radial and azimuthal d rift sp eeds in the 
test particle limit are (Nakagawa et ar][T986[ ): 



AVr 



AVa 



2t. 



1 



(57) 



1 



(58) 



where the dimensionless pressure support parameter, 

2 



V 



dP/dlnr 



1.5 



(59) 



and Vk is the Keplerian orbital speed. 

For loose coupling, azimuthal drift at the full AV^ « 
rjvia dominates the net motion. Eq. (55 1 will be satisffed 



(and the CTE will be ignorable) if turbulence is suffi- 
ciently strong so that SVg > rjVK- If we model the tur- 
bulence as SVg ~ \JaCg, then a> r\ ^ 10~^ is required, 
i.e. reasonably, but not exceptionally strong turbulence. 
For tight coupling, radial drift dominates relative mo- 
tion with AVr ~ V^kTs ^ rjVK- Eq. (551 is satisffed for 



SVg > rjVKT^ or a > rjr^. Thus even very weak tur- 
bulence {a > 10~^^ for mm-sized pebbles at 1 AU) will 
negate the CTE for Ts < 1. 

When the CTE is relevant, our results should apply 
approximately for an effective eddy time, t'^^^^, that is 
the smaller of the actual eddy time and the crossing time. 



t' 



eddy 



(V^cddy + l/^cross) 



-1/2 



*oddy/ \^ 



{AV/SVgr . 
(60) 

This correction factor is similar to that used by CDC in 
their equations (43) and (B12).^'' Our prediction also 
appears co nsistent with VJM R, in particular Fig. 1 of 
the related Volk et al. ( |1978 1 (where our relative drift 
AV is their VTT 



^® Vertical settling can be neglected as it is s maller than particle 
(and thus gas) RMS speeds from equation (|27|). 

Though CDC apply the correction for all St, even if tstop <SC 

eddy ' 

erated before crossing an eddy. 
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6.3. Particle Feedback and Concentration 

We model the passive response of particles to imposed 
turbulence, ignoring the aerodynamic feedback of par- 
ticles on turbulent gas dynamics. This restricts our 
analysis to local particle densities below the gas den- 
sity, pp <^ pg. This approximation is roughly valid for 
turbulence strong enough and/or particles small enough 
to prevent settling to such high densities. Specifically, 
Pp/Pg - i^p/^g){Hg/Hp) < 1 requires 

(61) 
En- 



10-2. 



from equation ( 23 1 , and a nominal Sp/Sg 
hancements of solid column densities relative to gas re- 
quires even stronger turbulence to remain in the test par- 
ticle limit. It is difficult to keep large particles from set- 
tling, but if only a small mass fraction resides in large 
particles, feedback effects might remain small. 

Drag feedback effects are known to trigger turbulence 
by themselves. Shear-induced instabilities from particle 
settling, i.e. Kelvin-Helmholtz instabilities, have been 
identified as a foil to the gravitational collapse of solids 
( Weidenschilling 1980[ CDC). Recently it has been shown 
that the in-plane drift of solids and gas triggers instabil- 
ities and t urbulence when feedback is included (Youdin 
and Goodman 2005; Yo udin and Johansen||200"7 ) . These 
"streaming" instabilities are particularly rele vant since 



they generate s ignificant particle clumping ( Johansen 
'and Youdin"2007) even when shear is also included (|jo- 



hansen et al. 2006). 

In addition to the streaming instability, other effects 
cause particles to concentrate in turbulent flows, as dis- 
cussed in the introduction. Most concentration mecha- 
nisms are size dependent and transient. (Appendix [B] 
mentions the likely destruction long-lived vortices in 3^ 
D). Nevertheless non-diffusive behavior, especially on the 
i ntegr al scale of turbulence as seen in [Johansen et al.\ 
( 2006 1 , may cause discrepancies with our simplified tur- 
bulence model. In the absence of gravitational collapse 
diffusion should still occurs on long length and time 
scales. Comparing our results to simulations (and pa- 
rameter regimes) with different degrees of particle con- 
centration will deepen our understanding of non-diffusive 
and other intermittent turbulent behaviors. For instance 
treating the global pressure gradients (and hence ra- 
dial drift) as a free parameter should modulate particle 
concentration. Drift that is either too strong (blasting 
through any eddy) or too slow (reducing the pileups in 
an eddy time) should reduce particle concentration. 

7. DISCUSSION 

We study the diffusion of particles, and the excitation 
of their random velocities, by turbulent forcing in Kep- 
lerian gas disks. We obtain an im prov ed re sult for ra- 
dial particle diffusion in equations (36 1 and (51), (iden- 
tical formulae from independent techniques)^ We ex- 
press radial particle diffusion in an ap proximate form 
for isotropic turbulence in equation \?>7\ and even more 
crudely (good for quick estimates) in equation (|5| , 



SCa: 



(1 



D 



p,2; 



1 -I- 4t2 



1 



(62) 



Our result improves on the estimate of particle diffusion 
in CDC, SccDC = 1 + St = 1-1- Ts/tc, by including or- 
bital dynamics and calculating diffusion for longer than 



a single eddy time (see J 4.1.1 1. The discrepancy with 
the CDC value is only significant for large particles with 
stopping times longer than the eddy or orbital time. 

O ur n ew results for in-plane velocity dispersions, equa- 
tion ( 33 1 , combine turbulent forcing and epicyclic mo- 
tions. Our results for vertical stirring in |4.2| agree with 
previous author s inclu ding CDC (see their eq. [53]) and 



DubruUe et al. (19951 for tight particle coupHng, and 
d'P for loose particle coupling. We add a correction 
for arbit rary eddy times which reduces particle veloci- 



ties (eq. 20 ) and scale heights (eq. 21 ) for long eddy 



times. We include simple, physically motivated prescrip- 
tions to generalize the scale height result to include thor- 
ough mixing of small particles with gas (eq. 28 ) and to 



describe the reduction of the effective (Lagrangian) eddy 
time due to particle drift (eq. [60]). 

We do not imagine that our results will significantly 
infiuence collisional velocities since the RMS velocities 
we calculate have a similar scaling to previous results. 
Order unity corrections are possible for collisions involv- 
ing marginally coupled particles, as orbital dynamics be- 
comes important and the drift through eddies is maxi- 
mized, but that calculation was beyond the scope of this 
work. 

The rates of gravitational collapse of large particles, 
which are limited by mass diffusion when Toomre's Q > 
1, will be increased. Preliminary investigations indicate 
that the rate of dissipative gravitational collapse remains 
constant as Tg ^ cx) while the CDC result would predict a 
growth rate that falls as (Youdin, in prep). Other in- 
stabilities limited by particle diffusion may be influenced 
as well. 

The detailed predictions of random velocities should 
prove useful in the interpretation of numerical simula- 
tions which couple particle dynamics to (usually MHD) 
turbulence. CFP already conflrmed their results for ver- 
tical oscillations with simulations of MRl driven tur- 
bulence. It remains to perform numerical tests of in- 
plane motions, but note that CFP have already recovered 
the expected epicyclic velocity dispersions (radial speeds 
twice azimuthal) for Tg = 10 particles. Ideally simula- 
tions of forced turbulence could investigate the effects of 
variable eddy times. 

We can apply our results the the chemical gradients in 
the solar system. A notable example is the spectroscopic 
zonation of the asteroid belt, but the issue is more general 
(Taylor 2001). We consider, very roughly, the minimum 



size a particle must have to drift radially faster than it 
diffuses, which is a crude measure of whether particles 
of a given size could maintain zonation over a drift time. 
Using the a-parameterization for tu rbu lent gas diffusion, 
Dg^x = o:xCg/f2, equations (l5| and (57) give the ratio of 



drift to diffusion time (over a scale of the disk radius 
r) as Adrift /idifi ^ a^/Ts. Any particle with > 
will diffuse slower than it migrates radially, i.e. a parti- 
cle bigger than a centimeter for w 0.01 at r ~ 2.5 
AU. This simple result deserves two comments. First it 
is the same criterion for a particle to be l arge enough 
to sediment to the midplane, from equation ( p3| (as long 
as ax ~ ctz which holds at least to order-of-magnitude) . 
Second, since ax ^ for subsonic turbulence, this fidu- 
cial particle size is tightly coupled, and could have been 
obtained equally well from the CDC result. Thus even 
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though CDC overpredict the diffusion of large particles, 
the global implications are minor: large particles drift 
aerodynamically faster than they diffuse (with either re- 
sult). 

The reduced diffusion of heavy particles would slow the 
loss of particles from pressure maxima (Whipple„1972j , 



including anti-cyclonic vortices, where radial drift is re- 
duced. Our improved result for diffusion is most relevant 
on smaller length scales, where diffusion rates are faster. 
This is the case for the dissipative gravitational instabil- 
ities mentioned above. 



APPENDIX 
PARTICLE SIZES AND DRAG REGIMES 



To translate from Tg to particle radius, a, we consider spherical particles of int ernal density p, — 2 g/cm^. We 
use a standard minimum ma ss nebula model ( Weidenschilling 1977 Hayashi 1981 1, which has a gas mean free path 
(evaluated in the midplane) ( Nakagawa et a?.||1986p , 



A 



Pgfmol 



1.2 m S 



(sAu) 



2.75 



(Al) 



where /j, = 3.9x10 g is the mean molecular weight, and (Tmoi = 2 x 10 cm^ is the molecular cross section ( Chapman 
and Cowling 1970). The gas surface density profile, Sg cx r^^^^ is normalized to U150 = ^g(5 AU)/(150 g/cm ) at 
r = 5 AU. 



When a < 9 A/4, Epstein drag from molecular collisions applies (Adachi et al. 19761. Then the particle size for a 
given (midplane) Tg is 

aEp(T,) = w 30 cm t,Si50 (^) ''^ • (A2) 



'27rp, 



For a > 9A/4, particles enter the viscous Stokes drag regime and 

ast(T-s) = 



9/iCgigtop 



4p,(Tinol 



90 cm ^/t. 



(sAu) 



5/8 



(A3) 



independent of gas density and with Cg as the sound speed. For linear drag the particle size is given by the smaller of 
OEp or ast- 

Drag forces become nonlinear when the Reynolds number. Re = AVa/{\Cg), of the flow around the particle exceeds 
unity. We take the relative motion of particles and gas, AV ~ rjv^ ^ Cg/uK, from the pressure supported gas rotation 
a small fraction of the Keplerian speed, = ^"f", see equation ( 59 1 . We estimate the size and stopping time, from 
equation (A3), of the transition to turbulent drag as: 



^s.turb 



aturb - AwK/cgW25 m S^r^^ (^) 

(sAu) 



4p.A 

9S^ 



2.5 



3.75 



! 700 S- 



(A4) 
(A5) 



Our assumption of a linear drag law applies up to large particle sizes in the outer disk, but only for Tg < 1 at r < 1 
AU. 

TURBULENCE MODEL FOR ARBITRARY EDDY TIME 



To numerically evaluate our results for in-plane particle stirring ([4.3 1, we relate the radial, azimuthal and correlated 
turbulent velocities ((Ug), (w^), and {ugVg}) as a function of Tq. The essential assumption is that velocities are isotropic, 
which we show is reasonable even for spatially sheared eddies with ^ 1. The detailed value of (ugWg) has a weak 
effect on particle stirring (and even then only near Tg = 1), but is included for completeness. We also argue that 
particle trapping is weak when eddies survive for only a single turnover time, even when that turnover time is long. 

We separately consider the limiting cases of large and small To. For To <C 1 we take the unsheared turbulence to be 
isotropic, (Ug) = (Wg). The velocity correlations are weak, and generated by the linear shear of radial velocity into 
azimuthal (which builds up for toddy) as (ugfg) ~ (3/2)To(Ug). 

The To 3> 1 case is more constrained by the dynamics of background shear. We introduce u' as the perturbed 
vorticity of the turbulence on the outer scale, which we compare to the background shear vorticity, wk = —3S7/2. 
(The ratio Icj'/tJKl is comparable to the ratio of turbulent speed to sound speed when the radial outer scale of the 
turbulence is comparable to the scale height.) We consider [w'/wkI < 1, because the efficiency of converting vorticity 
from the background shear into turbulence is likely to be less than 100%. (Note however that |w'/(jJk| > 1 is allowed 
for the forced To ^ 1 turbulence above.) 
To build a model for Ioj'/wkI < 1, we are guided by the dy namics of two-dimensional vortices, which have the 



following well-known properties (^Saffman [1992 Chavanis 2000): (1) vortices exist in regions where the perturbed 
vorticity lj' < is more negative than tEe background vorticity; (2) vortices are elongated in the azimuthal direction 
by a factor Ay /Ax relative to their radial extent; (3) the radial and azimuthal (with Kepler shear subtracted) 
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TABLE 2 
Symbols 



Symbol 



a 

^stop 
^cddy 
tsctt 

Ts 

Tc 

St 
SC(a;) 

"(2!,Z) 

Hp, _ffg 
SVg, 8ug, Swg 
SVp, 5up, Swp 

(UgVg),(UpVp) 



X, y, z 

Pp: Pg 

a,p, 

E, Eg 

V 
AV 



Definition/Use 



eqs. ( A2 I 



IA3I 



eq. Igi 
cq. 

■^^stop 
^stop /^ndH v 

gg/Dp, cq. |37| 
eqs. mj, ( [sip 
Dgn/c|7cq^^^ 
eq. |2l]| 

Tabrm 

eq. R FjLi] 



eqs. (10|l,l|Bll 

eqs. ( 13 
eqs. I Bl 




Meaning 



Keplerian orbital frequency 
particle stopping time 
eddy turnover time 
vertical settling time 
dimensionless stopping time 
dimensionless eddy time 

Stokes number 
Schmidt number (radial) 
diffusion coefficients: gas and particle (radial) 
dimensionless diffusion coefficients 
particle, gas scale heights 
RMS turbulent speeds: isotropic, radial, vertical 
RMS particle speeds: isotropic, radial, vertical 
mean squared turbulent gas speeds : 
radial, azimuthal, vertical 
mean squared random particle speeds : 
radial, azimuthal, vertical 
velocity correlations: gas, particle 
isothermal sound speed 
radial distance from star 
local radial, azimuthal, vertical coordinates 
particle, gas space density 
particle radius, internal density 
particle, gas surface density 
pressure support parameter 
particle-gas relative velocity 



fluid velocities are 



l^'^/wxly and Vg 



Thus at the vortex edge 



.2x1/2 



,2\l/2 



I Ax and the 



velocities have the same magnitude; (4) the time for fluid to circulate around a vortex is ~ 1/|^^'|- 

We now depart from our analogy with long-lived two dimensional vortices which are are absolutely stable in the 
absence of viscosity and are clearly not turbulence — since energy dissipation and angular moment um transport are 



absent. In three dimensions, though, vortices decay (Barranco and MarcusJ]2005| |Shen et al. 



20061. We hypothesize 



that three-dimensional turbulence is characterized by the properties (2) and (3) of vortices listed above; and that, based 
on property (4), vortices live for a time ~ l/l'^'l before decaying. Thus we assume teddy ~ l/l'^'l so the assumption 
jtj'/a^Kl ^ 1 is equivalent to Te ^ 1. 

To calculate the (ugUg) for Te ^ 1 we use the e stablished relation that the energy di ssipati on rate is 3l?/2 times 
the angular momentum flux (jLynden-Bell and Pringle 1974] [Lithwick and Chiang 20071. (See Youdin and Johansen 
2007 for how particle feedback introduces pressure work to the energy balance.) Since the energy dissipation rate is 
~ (itg -I- i'g)/toddy and the angular momentum flux is (ugVg), we have (ugVg) ~ 2(Mg -I- ■i;g)/(3re). 
Combining the results for the large and small regimes gives 

2 + Qr^ /4 

("g)-(^g)- M. (Bl) 

The numeri cal c oefficients should not be taken literally, but are for concreteness. The (ugVg) correlations given 
by equation (Bl I are a reasonable upper limit. In any event (ugVg) has only a moderate effect on particle stirring 

Sg, the standard approximation 

,2 



and only near Tg = 1. 



We also note that with an a-parametrization, (ugVg) = 
/aCg applies for hydro-turbulence only if correlations are strong (ugVg) 



I.e. near 



stresses, so Vt 



1 in equation (Bl I. In MRI turbulence on the other hand, a includes a dominant contribution from Maxwell 

\/aCg can hold with weaker correlations. 

Long-lived two-dimensional vortices can concentrate particles (de la Fuente Marcos and Barge 2001). However 



since we assume that eddies only survive for a single turnover time, particle concentration is weak, and our diffusive 
treatment should be a good approximation. Because the drift speed of a particle towards the center of a vortex is^^ 
Wdrift ~ VgTs/{l + Tg), the time for particles to reach the center of a vortex is Ax/wdrift ~ {^xhe){^ + Ts)/'''s 



Wdy(l + t^)/ts, i-e., it is always longer than icddy 
con centration to be importa nt . 



Vortices must survive for many circulation times for particle 



|Klahr and Henning 
is parallel to the midp 



19971 consider particle concentration by a "sideways" vortex, i.e., a vortex whose ro tation axis 



ane. Such a vortex concentrates particles in the particles' settling time if \uj\ < fi. Klahr and 



We neglect the extra drift caused by cause drag with the mean shear of the disk. This is appropriate if > H/r, because then 

the non-Keplerian velocities due to the vortex are larger than those due to the mean shear. If this inequality were violated, then particle 
would drift straight through vortices and concentration would be negligible. 
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Henning (19971 show that sideways vortices appear in 2D simulations of convection, so the efficiency of this mechanism 
in 3D witn radial drift deserves more study. 
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